Convergence of continuous-time quantum walks on the line 
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Abstract 



The position density of a "particle" performing a continuous-time quantum walk on the integer 
lattice, viewed on length scales inversely proportional to the time t, converges (as t tends to infinity) 
■ to a probability distribution that depends on the initial state of the particle. This convergence 

behavior has recently been demonstrated for the simplest continuous-time random walk ^ . In this 
brief report, we use a different technique to establish the same convergence for a very large class 
^>v^ I of continuous-time quantum walks, and we identify the limit distribution in the general case. 

I An article recently published in this journal 1 proves that a certain "continuous-time quantum 
walk" on the integer lattice Z has the same kind of convergence behavior as "coined quantum walks" 
t;;^- . on Z have ISIBl- In this note, the convergence result of P is generalized to a large class of continuous- 
I time quantum walks, using the techniques developed in for coined quantum walks. An interesting 
' aspect of the generalized continuous-time quantum walk is that the limiting measure need not have 
compact support. 



^ ■ Continuous time quantum walks on graphs were first defined in as follows. Consider a graph 

a ■ 

^ : 



^ , with vertex set V. Let denote the infinitesimal generator of the continuous-time Markov jump 
' process on V, where jumps only occur between adjacent vertices and the jumping rates are all equal 



to some 7 > 0. The continuous-time quantum walk of [Sj amounts to the unitary dynamics 



XI . V'o ^ e-**^-Vo (1) 

d I on £'^{V), the Hilbert space of square-summable complex- valued functions on V. Other authors [Z|IH1 
have defined continuous-time quantum walk as the dynamics 

^0 ^ e-**^Vo , (2) 

using the operator defined by the adjacency matrix A of the graph instead of L^. When the graph is 
the integer lattice Z, it hardly matters which way the continuous-time quantum walk is defined, for 
L.y = (1 — 2-y)I -\- 'jA and therefore and @ differ only by a change of time and phase. 

In m, Norio Konno studies the continuous-time quantum walk on Z defined as in © but with —t/2 
instead of t. That is, he studies the dynamics V'o e'^^^^'^^^ijjo, where A denotes the operator on ^^(Z) 
whose matrix with respect to the standard orthonormal basis {e„} is the adjacency matrix for the 
integer lattice (the standard basis vector is the member of £^(Z) with e„(n) = 1 and e„(A;) = for 
all k ^ n). In this case the operators e*^*/^^"^ can be expressed exactly in terms of Bessel functions, and 
the following convergence behavior becomes evident For each t, define the probability measures 

P,(n) = |(e„,e^(*/2)A^^u2 
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on Z. Then 



for -1 < a < 6 < 1. 



lim 

t— >cx> 



^u^u, Ja vrVl- 



at<k<bt 



This result is a special case of a much more general proposition. The main condition is that the 
generator A of the quantum walk Q be a self-adjoint operator that commutes with translations of 
£^(Z). Any such A has a matrix representation (with respect to the standard basis) of the form 
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(3) 



with ao = ao- We will only consider self-adjoint operators A that are Fourier transforms of a multi- 
plication operator, in the following sense. We denote the circle of unit radius by T and parameterize 
it by < 6* < 2tt. ^^(T, |f ) denotes the Hilbert space of square-integrable functions T. The Fourier 
transform 

(.F/)(n) = / me- 



-ine de_ 
2n 



is a unitary isomorphism from L^(T, ^) to with inverse 



We will assume that there exists a measurable real-valued function a{6) on T such that A satisfies 
{J^*AJ^f){e) = a{0)f{e) whenever / and af are both in L'^(T, We will further assume that a{9) 
is differentiable at almost every 9 € T. These assumptions on the form of A permit us vastly to 

generalize the main result of ^ fairly easily, but they are rather technical. Perhaps a more convenient 

00 

(but less general) condition is that the entries of the matrix Q satisfy Yl ^|an| < c«, for this implies 

n=l 

00 

that A has the desired form with a{6) = ao + X] (anC*"^ + a^e~*"^) continuously differentiable. 

n=l 

The main result of ^ is a special case of the following theorem: the case where a{6) = — cos{9) and 
V'o = eo- 

Theorem 1. Suppose that the matrix ^ represents an operator A on £^(Z) of the form TAi\a]T*, 
where Ai[a] denotes the operator on L^(T, ^) of multiplication by a measurable real-valued function 
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a{9). Suppose that a' (6) = ■^a{9) is defined almost everywhere with respect to Lebesgue measure dO. 
Let ipo be any unit vector in ^^(Z) and define 

Pt{n) = \{en,e-^'^i;o)f , (4) 
where Cn is the n*^ standard basis vector in £^(Z). Let Pi['i/'o] be the probability measure 

Ft[i^o]{dx) = Y,Pt{n)5{x-n/t) (5) 

on M, where 6{x — n/t) denotes the Dirac delta distribution at n/t. 

Then the probability measures Pt['i/'o] converge weakly as t — > oo to the probability measure PfV'o] 
defined on measurable subsets X cM. by 
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2tt 



F[i;o]{X) = — I , \{J^*i;omfde . (6) 



Proof: We will show that the "characteristic functions" of the probability measures ((5)) converge 
to the characteristic function of the probability measure for this implies the weak convergence 
of the probability measures themselves [HI Section XIII-1]. The characteristic functions ^t{^) of the 
probability measures Pt['0o] in Q are 

nez nez 
= (e-^^^V'o, E^/,e-''^i;o) = (V-o, e^'^E^/te-^'^i^o) , (7) 

where Ex denotes the multiplication operator {ExTp){n) = e*"^V(^) on £^(Z). To show that these 
characteristic functions converge, we will take Fourier transforms. 

Let A^[— a'] denote the multiplication operator {M.[—a']f){9) = —a'{9)f{6) and let 

H = J^M[-a]J^* . (8) 
Note that H may be an unbounded self-adjoint operator on £^(Z). We claim that 

lim e'^^E^ue-'^^i^ = e'^"^ (9) 

for all tf) € £^(Z). To prove this claim, first verify that 

{J^*e''^J^f){9) = {e''^*^^f){9) = e'''^^^^f{9) (10) 
{T*E^nTf){9) = f{9 + u/t) (11) 

for ah / G L'^{T, §). Supposing that {J^*ij){9) is a continuous function on T, we calculate that 

lim {J'*e'''^E^ue-''^iP){9) = lim {{J^*e''^r){J'* E^uJ'){T*e-''^r)J^*4^){9) 

lim e^*^(^)e-^*"('^+^/*)(^»(0 + uj/t) 



t— >oo 



-ium' (9) / T7* 



(12) 
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at almost every 9 by (|1U() and Since is bounded, the functions that converge pointwise in 

H12|) also converge in Lp'iT, ^) by Lebesgue's Bounded Convergence Theorem. The continuity of J- 
from L2(T, 1^) to implies that 

lim e**^-E^/ie-**^V = J^e^'^^[-^V> = e^'^'^V' 

as claimed in 0. This verifies the claim when is a continuous function; the general claim follows 
by a straightforward density argument. 

Applying the claim © in yields 

lim cD,(c^) = (^0, e^-^Vo) • 
Using the unitary isomorphism T and the definition ^ oi H one finds 

This is the characteristic function of the probability measure P[V'o] defined in (jUJ. □ 
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